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Abstract 

We consider a class of evolution equations in Lindblad form, which model the dy- 
namics of dissipative quantum mechanical systems with mean-field interaction. 
Particularly, this class includes the so-called Quantum Fokker-Planck-Poisson 
model. The existence and uniqueness of global-in-time, mass preserving solu- 
tions is proved, thus establishing the existence of a nonlinear conservative quan- 
tum dynamical semigroup. The mathematical difficulties stem from combining 
an unbounded Lindblad generator with the Hartree nonlinearity. 
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1 Introduction 

This paper is concerned with quantum mechanical multi-particle systems cou- 
pled to an external reservoir, i.e. so called open quantum systems |Dal BrPe . 
The dynamics of such systems can often be approximately described by ki- 
netic equations in the mean-field limit. Such self- consistent models appear in 
a wide range of physical applications, both quantum mechanical and classical, 
for example in gas dynamics, stellar dynamics, plasma physics, and electron 
transport. The corresponding nonlinear evolution equations are obtained as ap- 
proximations to the underlying (linear) many-particle models, and there exists a 
vast body of literature on their mathematically rigorous derivation: the classical 
Vlasov-Poisson system in [BrHel iBa) : the Hartree equation from the iV-body 
Schrodinger equation in the mean- field limit in |ErYa| : the Hartree-Fock equa- 
tion in |BaMa| . All of these models have in common that they fall into the class 
of Markovian approximation for the underlying dynamics and we refer to |Sp| 
for an extended overview of such derivations for a variety of kinetic equations. 
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In addition to a self-consistent Coulomb field we shall here be interested in 
quantum systems which in addition have a dissipative interaction with their 
environment. In many (practical) applications of such open quantum systems 
the interaction with a reservoir is described in a rather simple phenomenological 
manner, often using diffusion operators, quantum-BGK or relaxation-type terms 
|CaLellDeRillArT] when considered in a kinetic formalism. A prominent example 
of a linear open quantum system is the so called quantum optical master equation 
and its variants [GaZol IValj . However nonlinear mean- field models for open 
quantum systems also play an important role e.g. in laser physics (cf. |HeLi] 
and |Sp| for the Lieb-Hepp and the Dicke-Haken-Lax laser model, resp). 

In this work we shall be interested in a particular class of models which are 
frequently used in quantum optics DHR , iOCI IVal| and the simulation of nano- 
scale semiconductor devices |FMRI I.IuTaj , namely the quantum kinetic Wigner- 
Fokker- Planck equation (WFP) 

d t w + £,- \7 x w + e[V]w = Qw, x,^l J ,(>0, (1.1) 

which governs the time evolution of the Wigner function w{x,£,t) in (position- 
velocity) phase-space under the action of the potential V(x,t). In (|l.l(l the 
pseudo-differential operator Q[V] is defined by 

G[V}w(x,Z,t):=-^ [ I W(x+ V ~,t)-V(x- y -,tj 
(2ir) a J 7 K d xR d L V 2 J V 2 /J 

w(x,£,t) e ty - { Z-Z">d£'dy. (1.2) 

Q denotes the following diffusion operator 

Qw(x,(;):= D pp A s w + 2n div^w) + D qq A x w + 2D pq div x (V iW ), (1.3) 

with diffusion constants D (cf. (|2.17|l below) and the friction constant n > 
0. Here and in the sequel we set the physical constants H = m = e = 1, 
for simplicity. In semiconductor applications w(t, x, £) is the quasi-distribution 
of the electron gas and Q models (phenomcnologically) its interaction with a 
phonon bath. In our mean-field model the Hartree-type nonlinearity then stems 
from the repulsive Coulomb interaction between the electrons. Hence, <|l.lll is 
coupled to the Poisson equation 

AV = -n, (1.4) 

where n = J w d£ is the particle density of the electrons. 

Moreover, such Quantum Fokker-Planck (QFP) type equation are the most 
prominent model in the description of quantum Brownian motion, where a (mas- 
sive) quantum particle interacts with a heath bath and a possible external poten- 
tial, see e.g. [CaLel IDel IDil HaTl IQCj and |HuMa| . where this setting is proposed 
as a description of decoherence. Indeed most of these equations can be traced 
back to an early work by Feynman and Vernon |Fe Ve) . While formal deriva- 
tions of QFP equations were given in |CaLellDilllVaj . a rigorous derivation from 
many-body quantum mechanics is still missing, at least for the general class of 
models considered here. To the authors' knowledge, the only results in this di- 
rection are |CE FM FMRj, where special cases of the QFP equation arise, resp., 
in a space-time scaling limit and a weak coupling limit for a particle interacting 
with an infinite heat bath of harmonic oscillators, i.e. phonons. 
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In this paper we shall investigate well-posedness of QFP type equations with a 
mean-field Coulomb potential - the above mentioned Wigner- Poiss on- Fokker- 
Planck equation (WPFP) (|1 . 1|> — (|1 being one typical example. Specifically, 
we establish existence and uniqueness of global-in-time solutions to the Cauchy 
problem. Many of the analytical tools developed in the sequel will, however, 
directly apply to other open quantum systems in mean-field approximation (e.g. 
to the Dicke-Haken-Lax laser model) . First analytical results on the WFP and 
WPFP equations were obtained in SCDM (well-posedness of the linear 
equation, convergence to the unique steady state with an exponential rate), in 
|ALMS| (local-in-time solution for the mean- field model in 3D), and in |ACD] 
(global-in-time solution for the mean-field model in ID). 

In the mathematical analysis of mean-field QFP equations several parallel prob- 
lems have to be coped with: the Wigner framework often used in applications 
seems inappropriate since the particle density n = J w dt; is not naturally de- 
fined in this setup (typically, w G L 2 (R^xl^); cf. jATllALMSj for more details). 
We are hence led to study the equivalent evolution of the density matrix p(t) 
in the space of positive trace class operators J\. Moreover, in order to deal 
with the Hartree nonlinearity, an appropriate energy-space £ C J\ needs to be 
introduced, which is a generalization of the one used in BDF . In J\ the evolu- 
tion of the quantum system is then governed by a so called Markovian master 
equation, 



The considered Liouvillian C is obtained as a generalization of the one given 
by an inverse Wigner transformation of and will be stated in H2.5[) below. 
Since C (and in particular the included Lindblad operators [O]) are unbounded, 
this can be difficult even for linear equations and may lead to non unique and 
non conservative solutions. E.B. Davies showed in |Dalj that it is possible to 
construct, for a quite general class of unbounded Lindblad generators £, a so 
called minimal solution to the above master equation. However, this construc- 
tion is in general not unique, i.e. C does not uniquely determine a corresponding 
quantum dynamical semigroup (QDS) $>t(po) — e Ct Po- In particular, this im- 
plies that the minimal solution may not be conservative, i.e. trace preserving (cf. 
example 3.3 in |Dal| L which would be inappropriate for the above mentioned 
applications. 

While linear QDS have been studied intensively in the last three decades [FaRc 
I All I AlFa| . the literature on nonlinear QDS is no so abundant, see e.g. |Arll 
AlMe, BDF . By now, various sufficient conditions for the conservativity of 
linear QDS can be found in [ChFal |^GQ[ IHo| . For many concrete examples, 
however, these conditions are rather difficult to verify, as we shall discuss in more 
detail at the end of section |31 Moreover the assumptions on the nonlinearity 
introduced in |AlMej seem too strong for most physical applications. 

In this perspective, the present work establishes the existence and uniqueness 
of a conservative QDS for a concrete family of unbounded Lindblad generators 
L (including the WPFP model) with Hartree interaction. We shall consider 
Lindblad operators (representing the coupling to the reservoir) which arc linear 
combinations of the position and momentum operators, i.e. so called quasifree 
dynamical semigroups |Lilj . 




t > 0, 
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We briefly remark that the classical counterpart of WPFP, i.e. the Vlasov- 
Poiss on- Fokker- Planck system (and its linear version, the classical kinetic Fokker- 
Planck or Kramers equation [Rip 



d t f + ^V x f~V x V-V i f^DA i f + 2r 1 div^f), x,£eR d ,t>0 (1.5) 



allows for a much easier mathematical analysis. This is due to a natural 
L 1 (M^ x M^) -framework for 1)1.5(1 and to the positivity of the phase-space den- 
sity f(t,x,£), cf. |Bo| for the well-posedness analysis, |Ur| for existence of a 
unique steady state, and |DeVi| for convergence results to the steady state for 
the linear model. 

This paper is organized as follows: 

After introducing the model in section [21 we will prove in section [3] existence 
and uniqueness of a global, mass preserving solution to the linear equation, i.e. 
the existence of a conservative QDS. A crucial analytical tool towards this end 
is a new density lemma (relating minimal and maximal operator realizations) 
for Lindblad generators £ that are quadratic in the position and momentum 
operator. The mean field will then be included in section 0] (we shall restrict 
ourselves for simplicity to the case of d = 3 spatial dimensions). We prove that 
the self-consistent potential is a locally Lipschitz perturbation of the free evolu- 
tion in an appropriate "energy space" , and this yields a local-in-time existence 
and uniqueness result. Finally, we shall prove global existence of a conservative 
QDS in sectionElby establishing a-priori estimates for the mass and total energy 
of the system. 

2 The model equation 

In the sequel we shall use the following standard notations: 

Definition 2.1. J\ is the space of trace class operators on L 2 (M. d ) with the 
norm \\A\\ 1 := Tr where Tr denotes the usual operator trace on B(L 2 (H d )). 
J* C J\ denotes the subspace of self-adjoint trace class operators. Similarly, 

1/2 

Ji is the space of Hilbert- Schmidt operators with the norm \\A\\ 2 := (Tr \A\ 2 ) 
and II • loo denotes the operator norm in B(L 2 (M. d )). || • || 1 < p < oo is the 
norm of £P(R d )-functions. 

We consider open quantum systems of massive, spin-less particles within an 
effective single-particle approximation, as it has been derived for example in 
CEFM . Hence, at every time t E M a physically relevant, mixed state of our 
system is uniquely given by a positive operator p(t) S Sf(, in the sequel called 
density matrix operator. Since p is also Hilbert-Schmidt it can be represented 
by an integral operator p{t) : L 2 (R d ) -> L 2 (R d ), i.e. 



Its kernel p(-,-,t) E L 2 (M. 2d ) is then called the density matrix function of the 
state p and it satisfies lp(t)| 2 — I|p(vj*)ll2* By abuse of notation we shall 
identify from now on the operator p € 3{ with its kernel p(-, ■) E L 2 (M. 2d ). It is 
well known that we can decompose the kernel in the following form 




(2-1) 




(2.2) 



J'€N 
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where {Xj} & I (N) and the complete o.n.s. {V>j} C L 2 (R d ) are the eigenvalues 
and eigenfunctions of p. Using equation (|2.2I) one can define the particle density 
n[p] by setting x = y, to obtain 

n[p]{x) :=5>i |^-(x)| 2 , (2.3) 

However, since {x = y} C R 2d is a set of measure zero, this is not a mathemat- 
ically rigorous procedure for a kernel p{x,y) that is merely in L 2 (R 2d ). On the 
other hand, if p(x, y) is indeed the kernel of an operator p € J\ it is known, cf. 
|LiPa| . that the particle density can be rigorously defined by 

n yys P -M 2 /2e 



n[p](^):=Um/ p + j, x - ±) — — ^ dr, G L^(R d ). (2.4) 



And it satisfies = Tr(p) for p > 0. This issue of rigorously defining n[p] 

is one of the mathematical motivations for analyzing our mean field evolution 
equations as an abstract evolution problem for the operator p on the Banach 
space J( . 

Remark 2.2. Note that we can not use the decomposition 11'. 21 in order to 
pass to a PDE problem for the ipj, since the considered dissipative evolution 
equation in general does not conserve the occupation probabilities Xj. This is 
in sharp contrast to unitary dynamical maps generated by the von Neumann 
equation of standard quantum mechanics. 

We consider the following (nonlinear) dissipative equation modeling the motion 
of particles, interacting with each other and with their environment 

j t P = £(p):=-i[H,p]+A(p), <>0, 
p\ t=0 = po e Jl. 

Here, [•, •] is the commutator bracket, H and A(p) are formally self-adjoint and 
of Lindblad class. More precisely, we consider the Hamiltonian operator 

H:=~+V[p](x,t)-ip:[x,V} + , M e R, (2.6) 

denoting by [•, ■] , the anti-commutator. The operators x and V are, respectively, 
the multiplication and gradient operator on R d , i.e. [x, V], = x ■ V + V • x = 
2x ■ V + d. 

Remark 2.3. The operator H is sometimes called adjusted Hamiltonian, due 
to the appearence of the [x, V]+- term. Depending on the particular model, 
such a term may or may not be present, see e.g. |Del IDilj . Nevertheless it is 
included here, in order to keep our presentation as general as possible. 

The (real- valued) potential V is assumed to be of the form 

V[p}(x 7 t):=^ + V 1 (x)+ ( f>[p](x,t), xeR d , (2.7) 

where the first term of the r.h.s. denotes a possible confinement potential and 
V\ € L°°(R d ) is a bounded perturbation of it. We point out that the quadratic 
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confinement potential is not necessary for the subsequent mathematical analysis, 
it is just an option. <f> is the Hartree- or mean field-potential, obtained from the 
self- consistent coupling to the Poisson equation 

-A<f>[p]=n[p}. (2.8) 

For d = 3, we therefore get the usual Hartree-term: 

<t>[p}(x,t) = ±- ( x ,yeR 3 , (2.9) 

4tt J r3 \x - y\ 

where n is computed from p by (|2.4(l . This mean field approximation describes 
the (repulsive) Coulombian interaction of the particles with each other. 

The non-Hamiltonian part is defined as 

m 

A(p) := L JP L * 3 - 2 [ L i L i>p] + > m e N < ( 2 - 10 ) 

3=1 

or equivalently 

m 1 1 

^) = E 2 + 2 ^'"^ ' (2 ' n) 

i=i 

where the linear operators Lj (Lindblad operators) are assumed to be of the 
form 

Lj := a j -x + (3 j -V + i , a , fa € C d , j d e C. (2.12) 

Its adjoint is L* = 5>j ■ x — f3j ■ V + jj, and in the following we shall use the 
notation 

m 

/•: E 7 -: 7 - - ( 2 - 13 ) 

Remark 2.4. Linear models with Hamiltonians that are quadratic in the po- 
sition and momentum operator and with Lindblad operators of the form (|2.12(l 
give rise to so called quasifree QDS, and they are explicitly solvable in terms of 
Greens functions jLill ISCDM| . In oder to deal with nonlinear problems (in a 
"finite energy subspace" of J\) we shall, however, not use this representation, 
which moreover can not be generalized to higher order models, cf. remark [2. 71 

Remark 2.5. In the framework of second quantization and in d = 1, the space 
i 2 (R) is unitarily mapped onto ,F S (C), the symmetric or bosonic Fock space 
over C. This space is frequently used, for example in quantum optics, in order 
to describe two — level bosonic systems, cf. AlFa , GaZo . 
Assuming j = 0, /3 = 1 and a — 1/2, the Lindblad operators L, L*, become 
then the usual bosonic creation- and annihilation- operators 

af(x) := (| + d x )f(x), a* f{x) := (| - 8 x )f(x), (2.14) 

which, in contrast to the corresponding fermionic creation- and anihilation- 
operators, are unbounded. Of course, all results in our work can be equivalently 
interpeted in this framework of second quantization. 
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Example 2.6. A particularly interesting example in the above class is the 
Quantum Fokker-Planck equation (QFP). As a PDE for the kernel p(t,x,y) G 
L 2 (M. 2d ) it reads 



l9tp = \ 



-- + V(t,x),p 



P 



t=o 



+ iA(p), t > 0, 
p (x,y) G L 2 (R d x xR d ), 



(2.15) 



where 



A(fi) ■= - l(x - y) ■ (V a - V„) p + D M | V* + V y \ 2 p 

D pp . l2 21D M . . ,_ _ , (2-16) 

This model can be written in the form (|2.5(1 . (|2.1()(1 . iff the conditions 



2 

Z? ppJ D 9 , - ^ > !L ; ^ pp; > o, (2.17) 



hold (see |LilllALMS| for more details and a particular choice of the parameters 
p, ctj, f3j, jj). Using the Wigner transform |Wil [LiPa : 

W ( X ,U):=^l d p( X+ y -, X - y -,t)e^ d y. (2.18) 

the QFP equation (|2.15() can be transformed into the kinetic Wigner-Fokker- 
Planck equation (|l.ll) . In physical units D qq ,D pq ~ 0(H 2 ), cf. |Del IVa| . and 
hence we indeed obtain, at least formally, the kinetic Fokker-Planck equation 
(|1.5|l in the (semi-)classical limit K — * 0. Note that for 77 > 0, condition (|2.17|l 
implies that the diffusion operator Q from i|1.3fl is uniformly elliptic, which 
disqualifies the classical FP diffusion operator (i.e. D qq — D pq — 0) [KJ as an 
appropriate quantum mechanical equation. Nevertheless, this Caldeira-Leggett 
master equation |CaLe| is sometimes used in applications as a phenomenological 
quantum model, cf. |St j . 

Remark 2.7. To close this section we mention an interesting model from 
quantum optics which is not yet covered by our present analysis. The Jaynes- 
Cumming model with phase damping reads 

j t P = -i[H,p]+K[H[H,p}], (2.19) 

where k G M+ denotes the damping constant, cf. |Lo| . Since it involves Lindblad 
operators Lj that are quadratic polynomials of the position and momentum 
operators, it will be the focus of future research to (hopefully) extend the lemma 
13.71 (below) to such cases. 



3 Existence of a conservative QDS for the linear 
problem 

We consider the linear evolution problem on l 7 1 s (i 2 (K d )) 

Jt P = L ^ * >0 ' (3.1) 

p\ t=0 = po e Ji. 
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Here, C(p) := —i [H,p] + A{p) is the formal generator of a QDS on J(, with 

H = ~ + ^+V 1 (x)-i f i[x,V} + . (3.2) 

Definition 3.1. Given any Hilbert space TL, one defines a conservative quantum 
dynamical semigroup (QDS) as a one parameter Co - semigroup of bounded 
operators 

®t-Ji{H)^Ji{H), (3.3) 
which in addition satisfies: 

(a) The dual map $ t * : 13(H) -► B(H), defined by 

Tr(A<Z> t (p)) = Tr^*(A)p), (3.4) 

for all p €E J\(Ji\ A e B(H.), is completely positive. This means that the map 

$ t * ® I„ : ® B{H n ) <3 B(W„) (3.5) 

is positive (i.e. positivity preserving) for all n £ N. Here 7i„ denotes a finite 
dimensional Hilbert space and I„ is the n — dimensional unit matrix. 

(b) $4 is trace preserving, i.e. conservative (or unital). 

Remark 3.2. The notion QDS is sometimes reserved for the dual semigroup <f>j . 
Physically speaking, this corresponds to the Heisenberg picture. The appropriate 
continuity is then 

limTr(p($*(A)-A)) = 0, (3.6) 

for all p € Ji{Ti), A 6 B(H), i.e. ultraweak continuity. Complete positivity can 
be defined also for operators on general C*-Algebras A |Sti| and it is known 
that complete positivity and positivity are equivalent only if A is commutative. 
(Counter-examples can be found already for 2x2 complex valued matrices, see 
e.g. |AlFa| .) Again, from a physical point of view, complete positivity can be 
interpreted as preservation of positivity under entanglement. 



Following the classical work of Davies |Dal| we shall start to investigate the 
properties of the operator 

Y:=-iH-\l. (3.7) 

First we need the following technical lemma, the proof of which introduces some 
important notations used throughout this work. 

Lemma 3.3. Let P := pi(x, — iV) be a linear operator on L 2 (R d ) over the field 
C, where P2 is a complex valued, quadratic polynomial and specify its domain 
by 

V(P):={f:Ref,ImfeC%°(R d )}. (3.8) 
Then P is the maximal extension of P in the sense that 

V (P) = {/ G L 2 (R d ) : the distribution Pf e L 2 (R d )}. (3.9) 



Quantum Dynamical Semigroups 



!) 



Proof, (sketch) We define a mollifying delta sequence by 

<p n (x) := n d Lp(nx), x € K d , n € N, (3.10) 

with </? G Cg° and 90 > 0, = <p(— x), L d ip(x)dx = 1, supp<^ c {\x\ < 1}. 
Also, a sequence of radially symmetric cutoff function is defined by 

Xn{x):=x(^P), x€R d ,neN, (3.11) 

with Xn eC °°, < X < 1, suppx C [0, 1], x| [0f i] = 1- 

For / G L 2 (R d ) we define an approximating sequence in T>(P) by 

f n (x) :=Xn{x)(f*<p n ){x), neN. (3.12) 

We have to prove that for all / G L 2 (R d ), with Pf G L 2 (R d ), /„ -► / in the 
graph norm ||/|| p := ||/|| 2 + ||P/|| 2 . We clearly have 

/„"=5?/ini 2 (R d ), (3.13) 

and it remains to prove Pf n — * Pf in L 2 (R d ). This is now analogous to the 
proof of lemma 2.2 in |ACD| . when extended to complex valued functions /. A 
similar strategy is used again in the proof of lemma 1X71 below. □ 

Remark 3.4. Lemma 13.31 asserts that the minimal and maximal operators 
defined by the expression P = pt(x, — iV) coincide. This fact is closely related 
to the essential self-adjointncss of Schrodinger operators. The lemma provides 
an elementary proof of the well known fact that the Hamiltonian H = — A — \x\ 2 
is essentially self-adjoint on Cfi°(R d ), cf. corollary to theorem X.38 in |ReSi2| : - 

just apply the lemma to H with V(H) = C£°(R d ) and to H* ^ ^ . On the other 

hand, it is well known that H = —A + x 2 — x A is not essentially self-adjoint on 
Co°(R.), cf. example 1 of X.5 in RcSi2 . Therefore, lemma EOl can, in general, 
not be extended to higher order polynomials p(x, — zV). 

With the above lemma we can now prove that the main technical assumption 
on the operator Y (imposed in |Dalj . |ChFa| ) is fulfilled. 

Proposition 3.5. Let V± = and let the operator Y be defined on 

V{Y) := {/ G L 2 (R d ) : A/, \x\ 2 f G L 2 (R d )}. (3.14) 

(a) Then its closure Y is the infinitesimal generator of a Co - contraction semi- 
group on L 2 (R d ). 

(b) Further, the operators L jy L* : V{Y) -> L 2 (R d ) satisfy 

m 

(Yf,g) + (f,Yg)+Y / (L J f,L J g)=0, Vf,geV(Y), (3.15) 
i=i 

where (•, •) denotes the standard scalar product on L 2 (M. d ). 
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Proof. First note that for / E D(Y) the term x ■ V/, which appears in Yf, is 
also in L 2 (R d ). This can be obtained by an interpolation argument. Further, 
V(Y) is dense in i 2 (R d ), since C^(W l ) is. By Lemma IO we have 

V(Y) = {/ G L 2 (M d ) : Yf G i 2 (M d )}. 
Pari (aj: The proof proceeds in several steps: 

Step 1: We study the dissipativity of Y, which in our case is defined by 

Re(Y/,/)<0, V/6P(7). 
Since H from (|3.7fl is symmetric we obtain 

Re<iff/,/) = 0, V/gD(Y). 

Also we get 

-Re(L*L,/,f) = - (Ljf, Ljf) < 0, V/ e P(Y). 



Thus Y is dissipative and by theorem 1.4.5b of |Pa| also its closure Y is. 

Step 2: Its adjoint is Y* = iH — |i, with domain of definition T>(Y*). We have 
V(Y*) D V(Y), since 

(Yf,g) = {f,Y*g), Vf,geV(Y). 



As in step 1 we conclude that Y" 



is dissipative. We can now apply lemma 

V(Y) 



IO to P = F* with V(P) defined in El). Then P is dissipative on 

P(P) 

D(P) C D(Y) C D(Y"*). Since is closed, we have V(Y*) = V(P), the 
domain of the maximal extension. Thus Y* is dissipative on all of D(Y*). 

Step 3: Application of the Lumer-Phillips theorem (corollary 1.4.4 in jPaj ) to 
Y (with (Y) = Y*) implies the assertion. 

Part (b): We need to show: If/, Yf e L 2 (K d ), then Ljf, L*f G L 2 (R d ) follows. 
This can be easily seen from the fact that 



l^2{L j f,L j f)=-Be{Yf,f) < 



2 

3 



Equation H3.15[l is then obtained by a simple computation. □ 

With these properties of Y (as stated in proposition 13. 5fl . theorem 3.1 of Dal 
asserts that H3.lt has a so called minimal solution: 

Proposition 3.6. [Davies '77] There exists a positive Co - semigroup of con- 
tractions <I>t on . Its infinitesimal generator is the evolution operator C, 
defined on a sufficiently large domain T>{£), such that S\ 2 "D{£) 2 "D{7j). 
Here, Z : T>(Z) — > is the maximally extended operator with domain 

V(Z) = {pe J 1 s (L 2 (R d )) : Z{p) :=Yp + pY* G J? (L 2 (R d ))} . (3.16) 



From the above proposition we learn that the formal generator C, in general, 
does not unambiguously define a solution of the corresponding master equation, 
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in the sense of semigroups. Also, it is well known, that the obtained minimal 
solution need not be trace preserving (for nonconservative examples see e.g. 
|UaTllHop . 

On the other hand, if the semigroup corresponding to the minimal solution 
preserves the trace, it is the unique conservative QDS associated to the abstract 
evolution problem (|3.1|) , cf. |CGQ| IChFal IFaRel IHoj . We are going to prove 
now that in our case the minimal solution is indeed the unique QDS. To this 
end, we need to introduce some more notation: 

From now on we denote by 

(M(g)f)(x) := g(x)f(x), (C{g)f)(x) := (g * f)(x), g S C?(R d ), 

a family of multiplication and convolution operators on L 2 (R d ), where " * " is 
the usual convolution w.r.t. x. Further we define, for n € N, a family of sets 
V n C Jf{L 2 {R d )) by 

V n := {a n € J{ : 3p € J[ s.t. a n = M( Xn )C(cp n )p C(y>„)Af (*„)}, (3.17) 

where \n, 'fin are the cutoff resp. mollifying functions defined in the proof of 
lemma 13.31 above. For an operator p > with kernel (|2.2|l , the operator a n has 
an integral kernel given by 

<J n (x,y) = Xn(x)ip n (x)*p{x,y)*Lp n (y)Xn(y) 
x y 

= ^2*3 <Pj,ni x ) Vj,n(y)' ( 3 - 18 ) 
jen 

where <Pj, n (x) := Xn(x)(tp n * i/jj)(x) £ Qf(R d ) and ||<£j, n || 2 < = 1- Since 

a n > we get 

IIWIIi = Trcr n = ^\\fij,n\\l <J2^j = IIpIIIi- (3.19) 

jeN jew 

The union of all sets T> n will be denoted by 

2?oc := (J V n . (3.20) 

nGN 

Also we shall write for the graph norm corresponding to C 

\\ph~\ip\h + n£(p)lr ( 3 - 21 ) 

Then the following technical result, which is a key point in the existence and 
uniqueness analysis, holds. 

Lemma 3.7. Let V x = 0. Then: 

(a) The set is dense in J*. 

(b) C V(Z) C V(C) . 

(c) The operator C jj, is the maximal extension of C, in the sense that for 
each p € J", with C(p) € , there exists a sequence {tx n } n gN C 2?oo, such that 



lim \\p-(i n \\ c = 0. (3.22) 
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Proof. The proof is deferred to the appendix. □ 

Remark 3.8. For all p G 3f, £(p) can be defined (at least) as an operator 
C(p) : C^°(R d ) -> ^'( Rd ), the space of distributions. For £{p) G J? to hold, 
first of all an appropriate extension has to exist, such that £{p) G B(L 2 (R d )). 

We are now in the position to state our first main theorem: 

Theorem 3.9. Let V± = 0. The evolution operator £ generates on 3\ a conser- 
vative quantum dynamical semigroup of contractions $>t{p) — e ct p. This QDS 
yields the unique mild solution, in the sense of semigroups, for the abstract 
evolution problem 

Proof. Existence of $>t(p) — e p is guaranteed by proposition 13. 61 As a semi- 
group generator £ is closed, and by lemma 13.71 it is the maximally extended 
evolution operator. This implies uniqueness of the semigroup. Complete posi- 
tivity then follows from Stinespring's theorem jStil lAlFaj . 

It remains to prove the conservativity for the obtained QDS. This will be done 
by using a similar argument as in the proof of theorem 3.2 in Dal : 

Step 1: For the special case po G T>(£) the trajectory $>t(Po) is a classical solu- 
tion (in the sense of semigroups, cf. Eij), i.e. $>t(po) G C 1 ([Q,cx}),3i(L 2 (M d ))) 
and $ t (Po) € V(£), V t > 0. Hence Tr$ t (p ) G ^([0, oo),R) and we calculate 
for t > 0: 

^Tr$ t (p ) = Tr^M = Tr£($ t (p )) = 0. (3.23) 

To justify the last equality we note that Poo is || ■ \\c - dense in T>(£), by 
lemma T3.7I (c). Thus we can approximate $>t(po)i for every fixed t > 0, by 
an appropriate sequence {cr„} C Voo. Since T)^ is included in the domain of 
each "term" (IA.lt of the operator £ (as the proof of lemma ITT7T1 (b) shows), the 
cyclicity of the trace yields Tr£($>t(po)) = 0. Equation (|3.23|) then implies 

Tr$ t (p ) = Tr/5 =0, Vpo G V(£), t > 0. 

Step 2: The general case po G 3i{L 2 (M. d )) (i.e. <&t(po) is a mild solution) follows 
from step 1 and the fact that V(£) is dense in 3f{L 2 {R d )). □ 

From the above theorem, we obtain the the following corollary: 

Corollary 3.10. For p G T>(£) let 

£(p):=£(p)+£ p (p), (3.24) 

where 

£ p (p):=-i[V 1 ,p]+ £ L jP L*-±[L*L j ,p] + , (3.25) 

j—m-\-l 

with Vi G L°°(R d ), Lj E B(L 2 (R d )) and the sum converges in B{3i(L 2 (R d ))). 
Then the perturbed operator £ again uniquely defines a conservative QDS of 
contractions. 
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Proof. Existence and uniqueness of the Co-semigroup follows from standard 
perturbation results, cf. |Pa| . To prove conservativity of the perturbed QDS, 
let p(t) denote the solution of 

-j.P = £(p), p(o) = Po- 
The conservativity then follows from Duhamel's representation 

p(t) = $ t (p„) + / $ t _ s (£>(s))) da, (3.26) 
Jo 

by noting that Tr(£ p (p)) = 0. All other properties can be established by the 
same procedure as in theorem 1 of |AlMe| or by a Picard iteration. □ 

Remark 3.11. An alternative approach to prove theorem 13 . 91 could be to verify 
the sufficient conditions of ChFaj. In fact their assumptions Al and A2 are 
simple consequences of our lemma 13.31 and proposition 13.51 For their third 
condition A3 however, one would need to prove that Co°(K d ) is a core for Y 2 , 
defined on 

V(Y 2 ) := {/ G D(F) : Yf G V(Y)}. (3.27) 

With considerable more effort, the proof should be possible by extending the 
strategy of lemma l3~31 However, one can expect quite cumbersome calculations. 



4 Local-in-time existence of the mean field QDS 



We shall now prove existence and uniqueness of local-in-time solutions for the 
nonlinear evolution problem 



p = C{p), t>0 



d 

p(0) = P0 E J{. 
Here, the nonlinear map C is given by 
A 



+ V\p]-in[x,V] + ,p 



(4.1) 



(4.2) 



where the self-consistent potential V[p] is given as in 1)2. 7)1 and A(p) is the 
Lindblad operator defined by 1)2.10)1 and |2.12|) . 

To this end, we shall prove that the linear evolution problem 1)3.1(1 not only 
defines a Co-semigroup in J( (guaranteed by theorem 13. 9J) but also in an ap- 
propriate energy space. This is a parallel procedure (apart from severe technical 
difficulties) to solving the Schrddinger-Poisson equation in 7? 1 (R d ), cf. |GiVe| . 
Note that Davies' construction of a minimal QDS is valid only in J\. Hence, the 
required additional regularity of $t(po) has to be established explicitly. Also, 
one has to prove separately that this nonlinear model conserves the positivity 
and the trace of p. 

In the following, we shall restrict ourselves to the physical most important case 
of d = 3 spatial dimensions. 

Let us start by introducing the following definitions: 
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Definition 4.1. The kinetic energy of a density matrix operator p G is 
defined by 

E km [p] :=iTr(^=Ap^=A), (4.3) 

d 



where y— A denotes a pseudo-differential operator with symbol |£|, £ G R , i.e. 

^f^-T^df \£\Pf)(Z>* x <%, V/eff 1 ^). (4.4) 
Further, we define the external and the self-consistent potential energy of p € ^7/ 

by 

E-*[p]:=^Tr(Hp|x|), £ sc [p] := 1 TV(^]p). (4.5) 

The total energy will be denoted by 

E tot [p] := E km [p] + E ext [p] + E sc [p}. (4.6) 
In the sequel we shall work in the following energy space £ : 



£:={pe J? : Vl-A+MWl-A+M 2 G Jj 8 }, (4.7) 
equipped with the norm 



\\p\\ s := I yi-A + MWl-A+mii (4-8) 

This energy norm is a generalization of the one defined in |BDF| . In case p is 
indeed a physical state, i.e. p > 0, and if in addition p G 2?oo, one easily gets 

HpII £ = ||pIIi + IIV = ApV = a|Ii + IINpNIi, VpgPoo, p >o. (4.9) 

Hence, a density argument, similar to lemma 1X71 (c), implies for all p > that 
p G £ is equivalent to p € Jf and _E fcm [p] + £ e;E *[p] < oo. 

We further remark that in the above definitions we neglected the term —ip[x, V] + , 
which appears in the generalized (or adjusted) Hamiltonian operator H2.6[) of our 
system. Thus, even in the linear case, we have E tot [p] ^ Tr(Hp). The latter 
term would be the more common definition for the energy of the system. We 
note that we shall use E tot [p] only for deriving a-priori estimates and towards 
this end E tot [p] is the more convenient expression. 

Remark 4.2. Using the cyclicity of the trace, one formally obtains the more 
common expression for the kinetic energy of a physical state p > 0: 

E km [p] := ± Tr(V^ApV^A) = \ Tr(-Ap) > 0. (4.10) 

However, these two expressions for E km [p] are not fully equivalent, since Ap G 
31 requires more regularity on p than just requiring \f— Ap-\/— A G J*. (For 
more details see e.g. |Arj and the references given therein.) We further remark 
that if the kernel of p is given as in (|2.2I) the kinetic energy reads 



1 

j'ei 



£ fcin [p] = 2E A ^ v «^°- ( 4 - n ) 
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Similarly we get that for physical states p > it holds E ext [p] > 0, as well as 
E sc [p] > 0, since p > implies n[p] > and hence <f)[p] > 0, by l|2.9|) . Finally, 
note the additional factor 1/2 in front of the term E sc [p], which does not appear 
in the Hamiltonian Q2.fi [I . (|2.7J) . It is due to the self-consistent nonlinearity, cf. 

EEj. 

Using these definitions, we will now prove that the sum of kinetic and (external) 
potential energy is continuous in time during the linear evolution. 

Lemma 4.3. Let V\ = and po £ £, then 

(E km + E ext )[p{t)] S C([0,oo);E), (4.12) 

where p(t) := $t{po) G C([0, oo), J{) denotes the unique QDS for the linear 
evolution problem, given by yS.l)) . 

Proof. First, we note that each p £ £ C J( can be uniquely decomposed into: 
p = pi — p2, where 



Pi, 2 := A _1 (A ( oA) ± A _1 , A := ^/l-A + laf, (4.13) 

and (ApA^ denotes the positive resp. negative part of (ApA) £ J(. It holds: 
Pi, 2 > 0, as well as p\^ £ £■ Using this decomposition for the intial data po £ £ 
and since <i> t preserves positivity, we can restrict ourselves in the following to the 
case po > 0, hence p(t) > 0. The idea is now to derive a differential inequality 
for E km + E ext from ijO) . 

Let us define some energy functionals for positive p £ J{: 

E k k J[ P ] ■= --T1M), Hf\p\ ■= \ Tr(x k pxi), (4.14) 
with k, I = 1, . . . , d. For p £ X>oo, the cyclicity of the trace implies 

d d 

E km [p]=Y / E k k " k l [p], E^ t [p]=Y J Et%[p] (4.15) 
fe=i fe=i 

and, by a density argument, the formulas (|4.15|) also hold for p £ £. 

Step 1: We apply the operators x k , d k (from left and right) to (|3.1|1 and take 
traces. A lengthy but straightforward calculation, using the cyclicity of the 
trace and setting w.r.o.g. Trp(i) = 1, yields for the kinetic energy: 



d 

Tikin 

dt 2 

k=l k=l j=l fc=l k,l=l j=l 



d j 1 d m d dm 

E -oEE k>i 2 - v E - 2 E E Refo-k^iW 



- i ^ li'iui , ,vn ,,; ;. Tr(9 fc px ; ) + Im(a J:fc 7 i ) Tr(pd fc ) 
k ,i=i j=i 

+ i(^+J2 Tr (9kpx k ) S j ■ (4.16) 
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ex I 



For the external energy we obtain: 

d j 1 d m d dm 

E j t E k% = - 2 E E i 2 + 4 ^ E s m + 2 E E 

fc=l fc=l j=l k=l k,l=l j=l 

d m 

+ i E E M&.kiSw) + M/?i,fcT7) Tr(px fe ) 

fc,i=i j=i 

-i^+E^(S fc ^. (4.17) 

Step 2: These equations are not closed in E m and £ ,ext . To circumvent this 
problem, we shall use interpolation arguments: First, note that (dkpdk) G J\, 
iff {dkyfp) G J72, cf. |ReSil| . Thus we can estimate 

{Ml < IVpII IVpdkll = Mi IdkPdkh- 

Likewise, we get 

\dkPXi\\< l^fcVPll III VP^i III 2 = Wkpdkl-L llkwlli 

and one easily derives analogous estimates for the off-diagonal energy-terms 
^ext/kin ^ jj ence ^ es tiniating term- by-term in (|4.16|) . (|4.17|) . we finally obtain 



dt 

k=l 



k=l 



with some generic constant K > 0. Applying Gronwall's lemma then gives the 
desired result. □ 

This lemma directly leads to our next proposition: 

Proposition 4.4. Assume that po G £ and V\ G L°°(M. d ) s.t. additionally 
Wi G L«(R d ), for some 3 < q < oo. Then 

^t(po) G C([0,oo),£), (4.18) 

where $>t(po) denotes the unique linear QDS corresponding to 

Proof. The proof is based on a generalization of Grumm's theorem. As described 
in the proof of lemma FOI above, we only need to consider, w.r.o.g., the case 
p(t) > 0. 

Step 1: At first, one proves that for all /, g G L 2 (R rf ) and s > 0, 

hm(/,A P (i)A 5 ) = (f,Ap(s)Ag), (4.19) 

where (•, •) denotes the standard L 2 (R d ) scalar product. Choosing two sequences 
{fn}, {9n} C C °°(R d ), s.t. /„ ^ /, ffn "-^3° «? in L 2 (R rf ) the assertion then 
follows from a fairly standard approximation procedure. 

Step 2: Let Vi = first. By theorem 2.20 in jHIj (a generalization of Grumm's 
theorem), step 1 and the continuity of 

\jp(t)\t 1 + 2(E k ™+E™ t )[p(t)] = l\Ap(t)A\\ 1 
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(cf. lemma FO)) imply 

lim|||A (p(i)-p( s ))A|| 1= 0, Vs>0. 

t — >S 

Thus 

A P (i)Ae C([0,oo), J?(L 2 (R d ))) 

and the proposition is proved. The case V\ ^ can now be included by a 
standard perturbation result, cf. |Pa| under the additional assumption that 
VVi 6 L q (R d ), for some 3 < q < oo, cf. |Ar| for the detailed calculations. □ 

As a remaining preparatory step, the following lemma states an important prop- 
erty of the nonlinear mean field potential 4>\p\. 

Lemma 4.5. Let p £ £ and d — 3, then <f>[p] <G L°°(]R 3 ). Moreover, the operator 
[4>[p], p] is a local Lipschitz map from £ into itself. 

Proof. Once again we decompose p = p\ — p% s.t. pi^ > and pi^ € £, as 
given in (|4.13l) . In d = 3, we explicitly get from (|2.9[) 

1 x 

Therefore, the Hardy-Littlewood-Sobolev inequality and the generalized Young 
inequality, cf. |ReSi2| . imply for j — 1,2: 

4>\pj\ e Ll(R 3 )nL p (m 3 ), 3< P <oo, 

as well as 

V<p[pj] G L 3 J 2 {R 3 ) n L P {R 3 ), 3/2 < p < oo. 

Here, L p denotes the weak L p -spaces, cf. |ReSi2| . Hence, by a Sobolev imbed- 
ding, we obtain tf>[p] G L°°(R d ). Similar arguments as given in the proof of 
lemma 3.11 in |Arj then imply that [4>[p], p] is a local Lipshitz map in the energy 
space £ . To this end we first estimate 

HAWAII ! < IIA^A-iU llApAllx 

and use the assumption ApA € J\. For the first factor on the r.h.s. one calculates 
for / G Cg* 



2 



IIA^ma-VIII = II v (0MA- 1 /) 111- + \\VT+Wd>\p]A- l f\ 

We rewrite the operator of the first term on the r.h.s. as 

V MA- 1 ) = [((V0[ P ]) + 0[p]V)(l - A)- 1 / 2 ] [(1 - A) 1 /2 A - 1 

where both factors are in B(L 2 (R 3 )). The first factor is bounded since V0[p] G 
L 3 (R 3 ) and since (1 - A)" 1 / 2 is a bounded map from L 2 (R 3 ) into if^lR 3 ) ^ 
L 6 (R 3 ), due to a Sobolev imbedding. 

Summarizing we obtain 

\\m,p]\\ £ <C\\p\\l Vpe£, 
and the Lipshitz continuity then follows in a straightforward way. □ 
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We remark that the nonlinear map p i— > is continuous in £, but not in 

J( (L 2 (M 3 )) and this is the reason why we need to work in the energy space 
£ . However, the linear evolution problem (|3.1I) in general does not generate 
a contractive QDS on £ C Ji, except in the case of a unitary dynamic (i.e. 
Lj = 0). Hence, in order to obtain a global-in-time (nonlinear) existence and 
uniqueness result, we can not apply the results of |AlMej . which would require 
contractivity of the linear QDS in £. 

In the nonlinear evolution problem (|4.1|l the situation is even worse. Already 
in the case of a unitary time-evolution only E tot [p(t)] is conserved (for p — 
0), whereas ||/o(f)|| £ is not, due to the possible energy exchange between the 
potential and the kinetic parts. Hence a unitary but self-consistent evolution 
problem does not generate a contractive semigroup in £ either. 

With the above results, we are able to state the following local-in-time result: 

Theorem 4.6. Let p e £, d = 3 and V x £ L°°(R 3 ) s.t. Wi € L^(R 3 ), for 
some 3 < q < oo, then: 

(a) Locally in time, the nonlinear evolution problem j4-l\ ) has a unique mild 
solution &t(po) G C([0,T),£), where $*(•) denotes the nonlinear semigroup 
obtained by perturbing the linear QDS with the Hartree potential. This self- 
consistent potential satisfies: <j> € C([0, T); Ch(R 3 )). The map po i— ► &t(po) 
is Lipschitz continuous on some (small enough) ball {\\p — poWg < s} d £, 
uniformly for < t < T\ < T. Further, if the maximum time of existence T > 
is finite, we have 

Hm||d t (p )IU=oo. (4.20) 

t / i 

(b) For C(po) S £ we obtain a classical solution &t(po) £ C ([0,T),£). 

(c) The semigroup $ f is conservative. 

(d) The semigroup $t is positivity preserving and contractive on J7 1 S (L 2 (]R 3 )). 
Hence, it furnishes a nonlinear QDS: $ t : £ — * £ C J( . 

Proof. Part (a, b): By proposition 14.41 the unique conservative QDS $ t , ob- 
tained from theorem 13.91 also maps the energy space £ into itself. Lemma f4. 51 
and a standard perturbation result (cf. theorem 6.1.4 in |Pa| ) then yield the 
local-in-time existence of a solution for the nonlinear, i.e. mean field prob- 
lem. The continuity of <fi follows from the proof of lemma 1431 using $t(po) € 
C([0, T);£). The local Lipschitz continuity of the map po i— > ^t(po) follows 
from theorem 6.1.2 in |Pa| and the uniform lower bound for the existence time 
of trajectories $t(/o) that start in the neighborhood of po (cf. proof of theorem 
6.1.4 in |Paj). 

Part (c): The proof follows from Duhamel's representation, analogous to l|3.26|l . 

Part (d): Having in mind the result of part (a), we consider the nonlinear evo- 
lution problem l|4.1() as a linear evolution problem with time-dependent Hamil- 
tonian and write it in the following form: 

j t P = -i[H,p]+A(p)-i[<f>(t),p], t>0, 
p(0) = p > 0. 



Here, (j> € C([0, T); C(,(K 3 )) is the self-consistent potential <j>[p\. To prove the 
assertions of part (d), we shall approximate cf>(t) on [0,7i], T\ < T, by the 
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piecewise constant potential: 

0(t) := 0(t„), t n <t< t n+ i, < n < N - 1, 

with the uniform grid points: i n = nAi, At = T\/N, Hence, pit), t S [0, Ti] is 
approximated by Qv € C([0, Ti]; J? (L 2 (M 3 ))), solving 

-Qv = -i[ff,Qv]+^tev)-#(i),ftv], t>0, 
fi\r(0) = Po > 0. 

Since € Cb(IR 3 ), corollary 13.101 applies to the generator in (|4.22() on each 
time-intervall [t n ,t n -|_i]. In summary we have the following facts: 
4> is uniformly continuous on [0, Ti] w.r.t. || ■ H^, the solutions of 14.21J1 satisfies: 
|||p(i)||li < K, on < t < Ti, and the propagator corresponding to (|4.22|1 is 
contractive on J{{L 2 (]$?)). 

With these ingredients it is standard to verify that 

lim <; N = p, in C([0, Ti]; ^j*(i 2 (R 3 ))), 

TV— foo 

cf. the proof of theorem 1 in |AlMe| e.g. . Hence, the positivity of p(t) = $t(po) 
follows from the positivity of Sjsr(t). 

Analogously, the contractivity of the propagator corresponding to (|4.22|l implies 
the contractivity of $ t (po) in Jf(L 2 (R 3 )). □ 

Remark 4.7. If no confinement potential is present and Im^^fcOj,/) = 0, 
Vj, k, I, then theorem l4.6l also holds in the kinetic energy space £ km . In partic- 
ular, this is true for the QFP equation, where one can derive an exact ODE for 
the kinetic energy, cf. |A*LMS . 

In the next section we shall derive a- priori estimates on <&t(p) to prove the 
global-in-time existence of a conservative QDS for the mean field problem. 



5 A-priori estimates and global existence of the 
mean field QDS 

From theorem 14.61 we already know that ]|| j = |/OoIij f° r < < < T. It 
remains to prove an a-priori estimate on the energy of the nonlinear system. 
As a preliminary step, we introduce a generalized version of the Lieb-Thirring 
inequality: 

Lemma 5.1. Assume d = 3 and let p£ J{, p > be s.t. E kin [p] < oo. Then 
the following estimate holds: 

\\n\p]\\ P <K P IpI^M 1 - 9 , 1<p<3, (5.1) 

with 

<M> 

Proof. The proof is given in the appendix of |Arj , cf. also |LiPaj . □ 
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In the sequel this estimate will be used to derive an a-priori bound for the total 
energy. 

Proposition 5.2. Assume po S £, po > and d — 3. Then there exists a 
K > such that 

E tot [p(t)} < e Kt E tot [p ], < t < T, (5.3) 

where p(t) := &t(po), denotes the unique local-in-time solution of the nonlinear 
evolution problem 14-1\/ - 

Proof. Since $t is positivity preserving, we assume w.r.o.g. po > and hence 
have p(t) > 0, for all < t < T. The idea is again to derive a differential 
inequality for E tot . We first consider a classical solution <&t(po) € C 1 ([0,T),£) 
obtained from an initial condition with £(po) € £■ 

Step 1: We calculate the time derivative of the total energy, using the short 
notation p = j^p: 

= ^Tr (-iV=ApV=A + \\x\p\x\ +<Mp) ~ ^WiP) 
= Tr (-l^Ap^A + \\x\p\x\ + 4>[p]p^j + Tr(0[p]p) 

-~WWp). (5.4) 

For our classical solution p(i) the calculation (|5.4|l is rigorous since \\p\\e S 
C^O.T) and the self-consistent potential satisfies $ G C 1 ([0, T); C b (R 3 )). 

In order to simplify the last term on the r.h.s. of l|5.4|l we evaluate the trace in 
the eigenbasis of p (cf. I |2.2| l). This gives 

We now proceed as in |Arj : Integrating by parts several times and using the 
Poisson equation l|2.8|l . we obtain 

l±Tr{4>\p]p) = ~ [ \Vcj>[p](x)\ 2 dx = - I ^p}(x)A^[p}(x)dx 
i at i at j & 3 7 R 3 

j>[p](x)n{p](x)dx = Tr(</>[/>]/>). 

Inserting this into 1)5.4(1 . we get 

jE tot [p] = Tr Llv=Ap>/=A+ \\x\p\x\ + 4>[p]p 

= Tr Llv=aC(p)V=A+ \\x\C{p)\x\ 4 . ■[/>]/>) j . (45) 

In the following, we shall derive a differential inequality for E tot [p] from 
This expression is now considerable easier to deal with, since the self-consistent 
potential enters as if it was an additional external field (note that the factor 1/2 
in front of <p[p] has been eliminated). 
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Step 2: Similarly to the proof of lemma FOl we introduce an energy- functional 

Kt [P\ ■■= Ef? [p] + Elf [p] + [p] , k,l — 1, 2, 3, 
where E%f, E%? are defined as in (|4.14|) . Again, for all p G 2?oo, we have 

E tot [ P ]^jz E T k [p\ 

fc=l 

and, by a density argument, this carries over to p € £ ■ After some lengthy, but 
straightforward calculations (with extensive use of the cyclicity of the trace), 
we get from l|5.5|l . the following equation: 

3/33 \ 

i E^tfc = [j t E ~ \ E M^M)/' + (d k <t>[p])(d k p)) 

k=l V k=l fc=l / 

3 3 

+ J t E ^ + 2 ^E Tr(x fcP (a fc 0[p])) 

J fe=l fe=l 

3 m 

-i E ^2 lni ('^k( 3 j-j) ril ( x kP{di<P[p})) 

k,l=l j=l 
3 m 

- 1 E E Im (^,fc) Tr(p(9 fc 0[p])). (5.6) 
fe=i j=i 

Note that the first term of the r.h.s. of l|5.6[) - in big brackets - equals the 
time derivative of under the linear time-evolution. It is given by 14.16J1 . 
On the other hand, one easily checks that the time derivative of E^ x £ under 
the nonlinear time-evolution is equal to the linear one, hence given by 14.17|l . 
Since these kinetic and the external (potential) energy terms can be treated (by 
interpolation arguments) as in the proof of lemma 14.31 it remains to estimate 
the last three terms on the r.h.s. of Q5.6JI . 

Keep in mind, that we want to use a Gronwall lemma in the end. Hence, we 
need to find appropriate linear bounds for the r.h.s. of l|5.6Jl . (In the following 
we shall denote by K positive, not necessarily equal, constants.) 

Step 3: We first consider the term Tr(p(dk<f>[p]))- In order to calculate the 
trace, we need to guarantee that p{dk4>[p\) G Ji- Using the Sobolev inequality 
we estimate for ip £ L 2 (R 3 ): 



|(V=A + /)-Vlle < ^IKv^A + irVll^ < if IM 



2- 



since || (\/—A + I) ■ || 2 is an equivalent norm to || • ||^i. Holder's inequality and 
the bounds obtained in the proof of lemma l4~B1 then imply 

||(9^[ P ])(V^A + /)-Vll2<ll^[p]|| 3 ||(V=A + /)-Vll 8 

<K\\d k <f>[p]\\ 3 |M| 2 . 

In other words, (dk4>[p])(V— A + -0 -1 is a bounded operator on L 2 (K 3 ) and we 
get 

i\p(d k ct>[p])\h < ll(a fc 0W)(V^A + /)- 1 ILI (V=A + I)p\h 
<^||9^[p]|l3(^ fcin W + IIHi)- 
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Thus p{dk4>[p]) € J\i so we can calculate its trace in the eigenbasis of p and 
estimate it: 



IWW))I = 



<l|V^]|| 2 ||n[p]|| 2 . 



d k 4>[p](x)n[p\(x)dx 

The generalized Young inequality and the Lieb-Thirring inequality l|5.1|l imply 

\\Vm\2 < K\\n[p]\\ 6/5 < K\\pf^E kin \pfl\ (5.7) 
Further, using again Ij5.1|l . we have 

\\n\p]\\ 2 <K\\pi\ /4 E km [p]3/\ 
Hence, we obtain the following estimate: 

\Tv{p{d k ct>[p\))\<K\\p\\xE Un [ P ], (5.8) 

which is suitable for our purpose, due to the linear dependence on E km [p] . 
Step 4 : Next, we need to estimate the term 

3 

^2 6,* Tr ( x kp(di(l)[p])), 

k,l=l 

with the short-hand £k,i '■= ^ m (o^kPj,i)- To guarantee that Xkp{di<j)[p\) G J\, 
we only need to show y/p(di<p[p]) G J2, since we already know x^^fp G J2- 
This can be done as in step 3 above by noting that ^[p{\J — A + /) G J2 and 
(^A + Zl^WHeBfL 2 ^ 3 )). 

Hence, we can again calculate Ti(xkp{di(f)[p\)) in the eigenbasis of p: 

3 3 „ 

V £k,iTr{x k p(di4>[p])) = y2 tk.i x k di4>[p\(x)n[p](x)dx 

3 

= - Y] £k,i / a;fc9/0[p](a:)5^ </>[p](a:)cfe, 

(5.9) 



k,l m—l 



where we have used the Poisson equation l|2.8|l for the last equality. Integration 
by parts gives 

3 3 „ 

T Ck,iTr{x kP (d l tf>\p])) = £ 6,i / d l( f>lp}(x)d k <t>[ P ](x)dx 



k,l = l 



k.l=l 



3 

V 6,i / x k dl ' <j)\p]{x)d m <j)\p){x)dx. 

(5.10) 



k.l.m— 1 
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Adding the equations l|5.10|l and (|5.9|l yields, after another integration by parts: 

3 

2 Y, £k,lTr(x kP (di<t>\p])) 

k,l=l 

= E €k,l / di(f>d k (f>dx E £ k ,i / [x k d m (j)df m (j)-Xkdi4>dl im (j)]dx 

k,l=l 3 M ,m=l 3 

3 

= V] 6u / di(j>d k (j)dx 

3 

V / [<*fc,m <9? m <£ + ^ <^ m - 5ft,; <9^ m (j) - x k df <j>] <f>t 

3 

V &,* / |M 2 ^i ( 5 - n ) 



fc,Z,m— 1 
3 



= 2 E / 



di4>dk<t> dx 



k,l=l 



k,m—l 



where we write <f> = <f>[p] for simplicity and denote by Sk,i the Kronecker symbol. 
Therefore we can estimate 



E 6,iTr(x fc/9 (9^[p])) 



<#||V« 



where K depends on the coefficients £k,i- Hence, using the same estimates as in 
we have 



<K\\pll /2 E km [pY' 2 



<^Wi(Wlli + ^ m [p]), 

which is the desired linear bound. 

The third term in (|5.6() can be treated analogously to the previous case. 

Step 5: The steps 1-4, together with the estimates obtained in the proof of 
lemma FOl imply 

d 



j t E tot [p{t)] < KE tot [p{t)], 0<t<T, 



(5.12) 



with some generic constant K > 0. Applying Gronwall's lemma then proves the 
assertion. 

Strictly speaking, all the calculations of steps 2 — 5 first have to be done for an 
approximating sequence {a n } C 2?oo, such that a n — > p(t) in £ for each fixed 
t e [0,T) (cf. the proof of theorem l3.9fl . The estimate (|5.12|) then also holds for 
the limit p(t) since the constant K is independent of {<r n }. 

Step 6: So far we have proved l|5.3|) for classical solutions. By theorem |4.6f a'l any 
mild solution (i.e. <f>t(po) € C([0,T),£)) can be approximated in £ (uniformly 
on < t < T\ < T) by classical solutions. Hence <|5.3[1 carries over to all initial 
conditions P q S £ with po > 0. □ 

In view of (|4.20|l . and since ||/o(i)|| e < E tot [p(t)] we conclude from the above 
proposition that T = oo and obtain our main result: 
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Theorem 5.3. Let p £ S, d = 3 and V x G L°°(R 3 ) s.t. Wi G L ? (R 3 ), /or 
some 3 < g < oo: 

Then, the nonlinear evolution problem \4-l\l admits a unique mild solution, i.e. 
it generates a nonlinear conservative QDS: $t(po) € C([0,oo),£). 

6 Appendix: Proof of Lemma 13.71 

Without loss of generality we can assume that p is a nonnegative operator. (Oth- 
erwise one can split p into its positive and negative part |ReSil) and prove the 
result separately for each one.) Its eigenvalues are \j > and the eigenvectors 
are orthonormal. 

Part (a): For each p G J{ with finite rank N G N wc shall show that the 
approximation sequence {cr n } C 2?oo, defined in (|3.17() . satisfies a n — ► p in J\. 
With the kernel of a n G Jf* as in l|3.I8[) . we get from (|3.13() . that o~„ — ► p in the 
strong operator topology. Since we assumed that /? has finite rank, we conclude 
from l|3.19(l that the trace norms converge, i.e. 

lim \a n \ x = \p\ x . 

n — >oc 

Then Grumm's theorem (theorem 2.19 of [HI]) implies er„ — > p in J7i . 
Part f&J: The inclusion I?(Z) C T)(C) is already clear from proposition 13.61 
Thus it remains to show that for each a n G T> n c I?oo , with some fixed n G N, we 
have Z(a n ) G Jf*: First, note that Z(<7„) := ycrn + CnF* is a linear combination 
of the following terms (and their adjoints) 

x k a n xi, d k a n di, d k a n xi, x k xia ni d k di<j n , x k dia n , x k a ni d k a n , (A.l) 

where 1 < k, I < d and d k := d Xk . (Indeed not all of this terms really appear in 
the expression of Z ', but since the same argument for C is needed in the proof 
of theorem !3.9l we shall consider this more general case.) 

Since a n has a representation given by a n — M(x n )C(ip n )p C{ip n )M{x n ), for 
some p € S{, wc have to prove that the operator compositions x a V b M n C n are 
in B(L 2 (R d )). Here the multi-indices a, 6 G Nq are such that \a\ + |6| < 2. As 
an example we consider the operator x k di and write for / G L 2 (M d ): 

(x k diM n C n f)(x) = x k di{xn(x)(ip n * f)(x)) 

= X k [dlXn(x)((fi n * f)(x) + Xn(x){dl(fi n * f)(x)]. 

Since ip, x G C§° (see the proof of lemma rO| we have that 

|MM„C n /|| 2 < X M>n ||/|| 2 

and thus x k diM n C n G i3(L 2 (M d )). Hence x k d l a n = x k diM n C n pC n M n G Jf 5 . 
The other terms in (|A.1|I can then be handled in a similar way. 

Part (c): After the proof of part (a) it remains to show that for all p G J{ with 
C{p) G J{ , the following statement holds: 

lim lL{a n )-L{p)\ x =Q. 

n — >oc 

To simplify the proof, it is sufficient to consider a "model operator" JC(p), for 
which we choose I = k = 1 in (|A.1J| and further set all constants equal to 
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one. This simplification is possible since no cancellation occurs between the 
individual terms of K.(p). To simplify the notation further, we shall from now 
on write v := x\, d := d Xl . We choose K, in the form 

K{p) = Ki(p) + Ki{py, 

where 

/Ci (p) = vpv + dpd + dpv + v 2 p + d 2 p + vdp + vp + dp. 

The general (d - dimensional) case C(p) = —i[H,p] + A(p) described above is 
then a straightforward extension. The proof now follows again in several steps: 

Step 1: We write 

/CK) = K(M{ X n)C{y n )p C(<p n )M( X n)) 

= M{ Xn )C{<p n )K{p) C(<p n )M( X „) + Rn(p) + Rn{p)*- 

Since JC(p) <G , we can decompose it into JC{p) = JC + (p) — JC-(p), JC±(p) > 0. 
Applying part (a) of this lemma then yields 

lim I M{xn)C{y n )K(p) M( X n)C{y n )-K{p) ||, = 0. 

n^oo 

It remains to prove that R n {p) — » in J\, as n — > oo, which also implies 
R n {p)* — * in J\. For technical reasons (which will become clear in step 3) 
we split this remainder term into two parts: R n {p) = Rn(p) + Rn{p), an d treat 
each of them separately. 

Step 2: After some lengthy calculations, J?„(p) can be written as 

Rl(p) = M{dXn)C{iPn)pC{Vn)M{dXn) 

+ M(d 2 X n)C(<p n )pC(cp n )M( Xn ) + M(xn)C(vip n )pC(ip n )M(x n ) 
+ M(d X n)C(Lp n ) p C(<p n )M( Xn ) - M(x„)C(«Vn) P C(ip n )M( Xn ) 
-2M{x„)C(vLp n )pC( Vl p n )M(x„), 

where, on the level of the kernels, we have used several times the basic identity 
v (.f * 9) = v .f * 9 + f * v 9- Now we calculate for / e L 2 (M. d ) (remember v = x\) 

(C(xiip n )f)(x) := (x!- yi) ip„{x - y)f{y)dy 

JR d 

= I / n d+ \ Xl - Vl ) ip{n{x - y))f(y)dy = O (n" 1 ) . 

Thus we have |||C (vtp n )|| 00 = O (n^ 1 ) and similarly we obtain 

II|C(^)||| 00 = |||M( X „)||| 00 = 0(1), 

|||M(ax„)||| 0O = 0(n- 1 ), 

III^Vn)||| 0O = |||M(a 2 x n )||| 0O = 0(n- 2 ). 
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With these relations we can estimate 

III Rn(p) l 1 <|p| 1 |M(x„)Otf(«*'n)lL, 

+ MJM(x„)Otfte»)l JC'(''V, 

2 



+ IHIUIC^^I^IM^)!^ 

+ IIHIil|C(^)||Lll^(Xn)IL|M(9 2 x „)| 



l|p||llll|C(^)llloolll^(X»)ll| 2 ocll|C(^ 



il|C(^)||| 2 0O |M(x„)IL|M(9x r> 



= 0(n- 1 ) 



Thus Rn(p) — > uniformly in J\, as n — > oo. 

Step 3: Again a lengthy, but straightforward calculation shows that the second 
part of the remainder can be written in the form 

Rl (p) = M(ndxn)C( Vn )pC(^)M( Xn ) 

n 

+ M( X n)C(d(v<p n )) p C((p n )M( Xn ) + M(Xn)C(^) p C( Vn )M(nd Xn ) 

n 

+ M(nd X n)C(<p n ) p C(<p n )M(- Xn ) + M(xn)C(^) p C(nv<p n )M( X n) 

n n 

+ M(vdxn)C(ip n )pC(<p n )M(xn) + 2M(ndxn)C{^)pC(<p n )M( X n) 

n 

+ M(-Xn)C((p n ) pC(nv<p n )M(xn) + M( Xn )C(nvip n ) pC((p n )M(- X n) 
n n 

+ 2 M(- X n)C(nvtp n ) p C(<p n )M ( Xn ). 
n 

In contrast to step 2 these terms do not converge to zero uniformly in J\ , hence 
we shall proceed differently: 

As an example we consider the ninth term on the right hand side and write 

M( X n)C(nvtp n )p C(<p n )M(~ Xn ) = M{x n )C{nvip n )p N C{v n )M{- Xn ) 

n n 

+ M( X „)C(nv<p n )(p - p N )C(<p n )M(- X n), 

n 

where p N is the trace class operator p "cut" at finite rank N £ N, such that 
III 9 ~ P Hi — £ i £ e Direct calculations, similar to the one in step 2, imply 

I|C(n^ n )|L < K, |||M(n-^Xn)|IL < K, K £ K, (A.2) 
with K independent of n £ N. Thus we can estimate 

I|M( Xn ) C(nv<p n )(p - P N ) C^ n )M(-Xn)\h < sK\ (A.3) 

n 

Define II to be the projector on ran(p JV ). Then p N — Tlp N and 

\\C{nv^ n )p N \\ l < IC^^nlL \\p N \\ v (A.4) 

Now, since dim(ran(p JV )) < oo and since strong convergence equals uniform 
convergence on finite dimensional spaces [ReSilj , we get 

lim |||C(n U <p„)n||| 00 = 0. (A.5) 

n—>oo 
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Combining (|A.2|1 - (|A.5|) we thus have 

lim \\M(Xn) C(nv Vn )p N C(<p n )M(- Xn )h = 0- (A.6) 

n — >oo Yi 

Combining (|A.3Jl and l|A.6|) shows that 

\\M{ X n) C(nv^ n )p C(^„)M(-x„)| 1 

n 

can be made arbitrarily small for N sufficiently large. All other terms appearing 
in the expression of R\ can now be treated in the same way. 

In summary we have proved in steps 1 to 3 the assertion of the lemma. □ 
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